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Adsorption  of  Sulfate  on  Metal  Electrodes 


Leslie  A.  Barnes  Michael  R.  Philpott  and  Bowen  Liu 
IBM  Research  Diviaon,  Almaden  Research  Center 
650  Harry  Road,  San  Jose,  California  95120-6099 

Abstract 

The  adsorption  of  sulfate  (SO^')  and  bisulfate  (HSO4  )  on  copper  has  been 
studied  using  ab  initio  calculations  at  the  SCF  kvd  of  theory,  to  aid  in  the 
interpretation  of  in  situ  experimental  data  from  the  dectrochemical  interface, 
in  particular  optical  and  surface  X-ray  measurements.  The  calculations  are 
designed  to  ^ve  qualitative  insight  rather  than  quantitative  accnraqr.  Op- 
timixed  structures  and  harmonic  frequencies  and  inteudties  are  computed  for 
isolated  sulfate  and  bisulfate,  and  qualitative  agreemoit  with  experimental  data 
is  demonstrated.  The  effect  of  a  uniform  dectric  fidd  on  sulfate  is  also  stud¬ 
ied.  The  structures  of  and  Cu^HSOj  are  optimixed.  The  bonding 

is  dominantly  ionic,  with  consequential  small  orientational  dependence  of  the 
adsorbed  ligwd.  We  also  compute  structures  for  the  larger  CuJSOl''  duster. 
The  sulfate  ion  is  found  to  favour  a  “^own”  adsorption  geometry  with  sul¬ 
fate  oxygens  occupying  on-top  dies  over  the  previously  postulated  *‘l-down’’ 
structure  invoked  to  explain  the  early  surface  exafs  experiments. 


1  Introduction 

Understanding  the  adsorption  of  sulfate  on  surfaces  is  a  challenpng  sdentific  problem 
that  is  also  of  significant  tedmolo^cal  interest  in  connection  with  battery  chemistry, 
dissolution  of  metals  in  adds,  electroplating  and  polishing,  and  the  indpience  of  oxide 
films.  Recent  optical  and  surface  X-ray  measurements  of  metal  electrodes  in  con¬ 
tact  with  aqueous  sulfate  media  have  been  explained  with  hypothetical  models,  with 
sulfate  (S04'‘)  and  it*8  hydrolysis  product  bisulfate  (HSO4 )  in  spedfic.  orientations 
adsorbed  on  the  metal  surface  (IH6]-  Figure  1  shows  some  of  these  modds. 

The  high  symmetry  (Tii)  and  many  vibrational  modes  of  sulfate  that  are  exdu- 
sivdy  IR  or  Raman  allowed  make  it  a  potentially  useful  probe  of  forces  influendng 
orientation  and  of  the  local  field  within  the  electrochemical  double  layer.  Bisulfate, 
with  lower  symmetry,  has  more  vibrational  frequendes  which  may  be  IR  or  Raman 
allowed,  some  of  which  overlap  with  sulfate  modes.  This  jomplicates  interpretation 

’Mailing  Address:  Mail  Stop  BTC-230-S,  NASA  Ames  Research  Center,  Moffett  Field,  California 
94035-1000 
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b.  thre«  down 


Figure  1:  Proposed  Sexafs  model  of  sulfate  on  upd  Cu  monolayer  on  Au(lll) 
surface  along  with  some  optical  models  for  bisulfate. 

of  the  experimental  data,  where  both  sulfate  and  bisulfate  may  be  present.  ¥tom  a 
computational  viewpoint,  the  larger  size  of  sulfate  and  bisulfate  make  the  calculation 
and  interpretation  of  structures  and  spectra  more  difficult  than  for  smaller  ligands 
sudi  as  carlxm  monoxide  or  cyanide  [7].  Nevertheless,  because  of  the  existence  of  sev¬ 
eral  piv^KMed  structxires  on  Cu  surfaces  (5, 6]  and  the  variety  of  other  data  bearing  on 
adsorption  site  geometry  and  structure  (SHU]*  if  i*  clear  that  model  ab  initio  quan¬ 
tum  diemical  calculations  could  play  a  crudal  role  in  duddating  gecmietiy,  structure 
and  bonding  at  the  adsorption  site. 

The  goal  of  the  present  work  is  to  use  a  series  of  simple  modds  to  explore  the 
interaction  of  sulfate  and  bisulfate  vdth  metd  surfaces  and  aid  in  the  interpretation 
of  tn  min  ejqperiments.  Previous  calculations  have  studied  Cu  dusters  with  smdl 
ligands  like  0,  CO,  and  CN"  On  ■ilver  dusters  the  ligands  F*",  Cl~,  Br“, 

I~,  azide  Nj,  and  thiocyanate  SCN‘  have  been  studied  (idHlS].  The  polyatomiv 
and  potentially  polydentate  ligands  S04~  and  HSO4  are  more  demanding  of  compu¬ 
tational  resources  due  to  the  larger  number  of  atoms  and  dectrons  and  the  variety 
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performed  [10].  Other  structure  determinatioos  have  used  STM  [30]  and  AFM  [31] 
to  obtain  geometric  information  about  the  effect  of  sulfate  on  metal  adatom  sites  on 
gold  surfaces. 

Optical  spectroscopy  has  provided  important  information  about  adsorption  of 
bisiilfate  and  sulfate  on  metals  and  also  how  vibrational  frequencies  change  on  ad> 
sorption  and  with  changes  in  electrode  potential  In  particular  in  sUu  surface  FTIR 
spectroscopy  [lH^]«  shown  how  sulfate,  bisulfate  and  water  adsorb  on  poly- 

crystalline  and  single  crystal  surfaces  and  single  crystals  with  submonolayer  of  under 
potentially  deposited  metals.  Finally,  we  point  out  that  other  optical  spectroscopies 
sudi  as  second  harmoric  generation  and  sum  frequency  generation  [33]  have  provided 
some  information  concerning  the  interaction  of  adsorbed  sulfate  on  metal  electrodes 
over  a  wide  range  of  potentials  [34,  35]. 


3  Methods 

For  isolated  SO4'  and  HSO4  we  cany  out  ab  initio  all-electron  SCF  calculations 
using  analytical  derivative  methods  to  determine  optimal  geometries,  harmonic  fre¬ 
quencies,  and  both  IR  and  Raman  intensities.  The  SCF  method  is  expected  to  give 
bond  stretches  which  may  be  10-20%  too  high  compared  to  experiment.  However,  the 
results  should  be  qualitatively  correct.  For  SO4”,  we  have  also  carried  out  geometry 
optimizations  and  frequency  calculations  in  the  presence  of  an  applied  electric  held. 
These  w^  carried  out  under  the  constraint  that  the  S  position  was  fixed.  We  have 
also  computed  IR  intensities  for  S04~  in  the  presence  of  an  applied  field,  but  not 
Raman  intensities.  For  the  single  copper  atom  and  copper  cluster  calculations  with 
sulfate  or  bisulfate,  we  carry  out  ab  initio  all-electron  SCF  calculations  using  analyt¬ 
ical  derivative  methods  to  determine  optimal  geometries,  with  the  constraint  of  Csv 
symmetry  for  Cu^SOj”  and  C,  symmetry  for  Cu^HSO^.  The  CujSOj”  calculations 
were  constrained  so  that  the  Coi  geometry  was  fixed  at  that  of  bulk  Cu  metal  on  the 
(111)  surface  [36]  and  the  CujSO*"  cluster  had  C3,  symmetry. 

The  sulfur  basis  is  derived  from  the  (12s  7p)  primitive  Gaussian  set  of  Dunning 
and  Hay  [37],  contracted  to  [6s  4p]  and  supplemented  vdth  two  d  functions  with 
exponent  0.9  and  0.3.  The  oxygen  bi^is  is  the  (9s  5p)  primitive  Gausrian  basis  set  of 
Huzinaga  [38],  contracted  to  [4s  2p]  according  to  Dunning  [39].  This  is  supplemented 
with  a  diffuse  p  function  vrith  exponent, 0.059  [37]  and  a  d  function  with  exponent 
0.72  [40].  The  basis  for  hydrogen  is  the  (4s)  primitive  Gaussian  set  of  Huzinaga, 
contracbMl  to  [2s]  according  to  Dunning  and  Hay,  and  supplemented  with  a  single 
p  function  with  exponent  1.0.  For  Gu  we  use  the  (14s  9p  5d)  primitive  Gaussian 
basis  set  of  Wachters  [41],  contracted  to  [8s  4p  3d]  using  his  contraction  scheme  2. 
Two  diffuse  p  functions,  as  recommended  by  Wachters,  and  the  diffuse  d  function  of 
Hay  [42]  are  added,  yielding  a  final  basis  set  of  the  form  (14s  lip  6d)/[8s  6p  4d}.  Inis 
basis  should  be  adequate  to  describe  the  dominantly  electrostatic  interaction  between 
the  Cu  4s  electron  and  SO*”  or  HSO7  units. 
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I  of  potential  binding  modes.  We  first  examine  the  isolated  sulfate  and  bisulfate  ions, 

and  then  look  at  the  effect  of  a  strong  electric  field  on  sulfate.  We  then  model  the 
binding  to  the  surface  first  with  a  single  Cu  atom,  and  then  progress  to  larger  clus¬ 
ters  with  four  metal  atoms,  looking  at  geometries  and  relative  binding  ener^es,  using 
all  electron  ab  initio  SCF  Hartree-Fock  calculations  and  analytical  derivative  meth¬ 
ods.  A  more  complete  description  of  additional  calculations  for  sulfate  and  bisulfate 
including  adsorption  on  ten  atom  Cu  clusters  will  be  reported  elsewhere  [19]. 

This  paper  is  organized  in  the  following  way.  In  §  2  we  provide  a  survey  of  the 
experimental  data.  In  §  3  we  discuss  the  theoretical  methods  used,  and  in  §  4  we 
present  our  results  and  discussion,  first  presenting  results  for  S04~  and  HSO^,  then 
Cu‘’SO^,  Cu®HS07  and  Cu$SO^. 

2  Survey  of  Experimental  Data 

In  this  section  we  briefly  describe  some  of  the  more  important  experimental  work  that 
this  paper  impacts.  There  are  three  main  areas:  adsorption,  structure,  and  optical 
spectroscopy. 

In  aqueous  solution  sulfate  belongs  to  the  class  of  anion  known  to  undergo  specific 
adsorption  [20,  21]  in  the  range  of  electrode  potential  corresponding  to  a  Uietmody* 
namically  stable  dectrodhenoical  double  layo*.  In  this  process  more  weakly  bound 
water  molecules  between  the  electrode  and  the  ion  are  displaced  and  the  ion  makes 
physical  contact  with  the  surface  atoms  of  the  metal.  This  is  the  origin  of  the  term 
^contact  adsorption”.  Early  discussions  of  this  phenomena  also  considered  whether 
a  chenucal  bond  was  formed,  but  conventional  wisdom  now  describes  the  adsorption 
as  primarily  physisorption.  In  practice  contact  adsorption  results  lErom  a  balance  of 
opposing  free  energy  interactions  of  the  following  type:  dectrostatics,  enthalpies  of 
hydration  of  ion  and  surface  atoms,  and  entropy  of  displaced  water  molecules.  G«i- 
erally  speaking,  large  ions  (eg.  I”,  Cs^)  contact  adsorb  whereas  small  ions  (eg.  F~, 
Ca'*')  do  not.  The  contact  adsorption  series  for  noble  and  coinage  metal  electrodes 
immersed  in  aqueous  electrolyte  is  BFj’  <  PF?  <  ClOlf  <  F”  <  HSO4  <  SO*"  < 
Cl"  <  SCN"  <  Br"  <  I"  [21,  22).  Verification  of  contact  adsorption  for  sulfate  and 
bisulfate  on  metals  has  been  demonstrated  by  a  variety  of  classical  electrochemical 
techniques  including  cyclic  voltammetry  [23H25],  differentia]  csq)adtance  [26],  and 
radio  tracer  analysis  [27]-[29]. 

Structure  detemunations  by  surface  X-ray  scattering  have  proved  very  important 
because  actual  geometric  information  has  been  derived.  For  example  Blum  et  al.  [5, 6] 
proposed  the  first  structure  for  sulfate  adsorption.  These  results  are  particularly 
important  because  it  was  proposed  that  sulfate  may  occupy  an  on-top  site  with  the 
Gu-O-S  bond  parallel  to  the  normal  of  the  (Ill)  surface  plane.  Oomplimeutaiy  X- 
ray  studies  have  been  reported  for  various  Au  surfaces  [8,  9].  Similar  on-top  site 
geometries  have  been  proposed  for  adsorption  on  Cu(lOO),  obtained  by  emersing  the 
electrode  from  acidic  sulfate  solution  prior  to  tranfer  into  UHV  where  LEED  was 
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Figttre  2:  Isolated  S04~  and  HSOT  structures.  Bond  lengths  are  in  a.u.  and  angles 
are  in  degrees  (1  a.u.  s=  0.529177  A). 

The  programs  CADPAC  [43]  and  the  Cambridge  Direct  SCP  [44]  axe  used,  im* 
plemented  on  an  IBM  3090/300J  and  RISC  SYSTEM/6000  computers  at  the  IBM 
Almaden  Researdi  Center.  CADPAC  was  modified  to  compute  frequencies  and  inten¬ 
sities  in  the  presence  of  an  applied  electric  field  according  to  the  methods  of  Duran  tt 
<d.  [45]. 

4  Results  and  Discussion 

4.1  Isolated  SOj'  and  HSOJ  ions 

To  illustrate  the  qualitative  accura/^  of  the  calculations,  we  consider  here  the  struc¬ 
ture  and  vibrational  fipequendes  of  isolated  804“  and  HSO7.  We  emphaaze,  however, 
that  our  results  are  for  ^gas-phase”  spedes  —  so  the  comparison  with  mqperiment  will 
only  be  qualitative.  Neverthdess,  the  results  axe  useful. 

The  structure  of  isolated  S04~  and  HSO7  bk  illustrated  in  Figure  2.  Isolated 
S04~  has  tetrahedral  symmetry  (TV  point  group),  whereas  HSO7  has  a  single  plane 
of  symmetry  (C«  point  group).  The  computed  ^nd  distance  for  S04~  of  2.80  a.u.  is 
in  good  agreement  with  the  experimental  value  of  2.82  a.u.,  which  is  the  value  in  the 
Na3S04  crystal  and  also  is  the  average  value  in  the  K3SO4  crystal  [36].  The  computed 
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Table  1:  Siunmaxy  of  the  SOj"  harmonic  frequendes  and  Intensities 


Mode 

FVeq. 

cm“^ 

IR  int. 
km/mol 

Raman  int. 
AVamu 

Expt.  FVeq.“ 
cm“* 

471 

0 

2.3 

451 

666 

50 

2.4 

618 

1035 

0 

24.6 

981 

1174 

605 

8.2 

1104 

•  Aqueous  Na2S04  [49] 


structure  of  HSOJ  is  also  in  general  agreement  with  data  &om  crystal  structures  of 
KHSO4  [46]  and  NH4HSO4  [47].  For  example,  in  KHSO4  the  S-03  distance  is  around 
2.72  a.u.,  whereas  the  S-01  distance  is  longer  at  around  2.77  a.u.,  due  to  hydrogen 
bonding  with  other  HSOJ  units  in  the  crystal.  The  S-02  distance  is  around  2.96  a.u. 
in  the  crystal,  shorter  than  the  theoretical  result,  again  an  effect  due  to  hydrogen 
bonding  in  the  crystal.  The  0-H  distance  is  only  1.37  a.u.  in  KHSO4  and  NH4ELS04, 
compared  with  the  1.79  a.u.  value  found  here.  This  difference  is  larger  than  could 
reasonably  be  expected.  However,  as  discussed  by  Nelmes  [47],  the  “expected”  0-H 
bond  distance  in  other  hydrogen  bonded  systems  is  in  the  range  1.8~1.9  &.u.  For  ex¬ 
ample,  the  0-H  distance  in  H3SO4  [2, 48]  is  1.83  a.u.  Therefore  the  theoretical  value 
seems  very  reasonable.  Hoarding  the  an^es,  there  is  again  qualitative  agreement, 
except  for  the  ZS-O-H  which  is  not  unexpected  based  on  the  bond  distance  differ¬ 
ences.  For  example,  in  KHSO4,  the  Z03-S-02  is  around  105”,  compared  with  104“ 
theoretically,  and  the  Z03-S-01  is  111*  compared  with  115*  theoretically.  There  is 
some  difficulty  making  direct  oompansons,  due  to  hydrogen  bonding  and  other  “crys¬ 
tal”  effects  —  however,  the  theoretical  results  are  qualitatively  correct.  Thus,  for 
HSO7,  the  lengthening  of  the  S-OH  bond  distance  and  the  shortening  of  the  other 
S-O  distances  compared  ^th  the  S-O  cUstance  in  SO’**  is  seen  in  both  the  exper¬ 
imental  and  theoretical  results,  and  the  changes  in  the  bond  angles  are  also  in  the 
same  general  direction. 

The  computed  frequencies  and  intennties  for  SO4'  are  ^ven  in  Tsbie  1  and 
are  in  good  agreement  with  the  aqueous  Na2S04  results,  although  the  symmetric 
stretch  is  too  hi^,  as  expected  at  the  SCF  kvri  of  theory.  In  general,  firequendes 
which  are  too  high  also  cmreiate  with  bond  distances  which  are  too  short  at  the 
SCF  levd  cff  theory  —  however,  external  media  effects  in  the  solution  or  crystal 
environment  make  spedfic  comparison  between  the  hequendes  and  bond  lengths 
difficult.  In  fact  the  agreement  between  the  computed  and  experimental  frequendes 
is  somewhat  fortuitous  —  the  aqueous  results  will  undoubtedly  be  affected  by  factors 
such  as  hydrogen  bonding  [49],  whereas  the  effect  of  electron  correlation  has  not  been 
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Table  2:  Comparison  of  computed  HSO4  frequencies  and  int^isities  with  experiment 


■ 

Ereq, 

cm'* 

KiSiBI 

a" 

93 

2.3 

145 

a' 

441 

4 

1.5 

447 

a" 

476 

11 

1.1 

418‘ 

a' 

616 

52 

2.8 

576 

a" 

641 

71 

2.8 

589 

a' 

642 

27 

3.8 

610 

a' 

861 

330 

9.5 

883 

a' 

1157 

123 

18.4 

1022 

sf 

1258 

132 

3.0 

1044 

a" 

1351 

521 

4.7 

1194 

a' 

1409 

467 

4.0 

C 

a^ 

4176 

73 

61.2 

3510 

*  FVequendes  from  Raman  spectrum  of  NH4HSO4  crystal  (See  reference  [49]  and 
references  therdn) 

*  There  is  another  a"  mode  reported  at  407  cm“* 

‘  FVequendes  in  the  range  ISOO-lSOOcm"^  have  been  reported  from  IR  spectra  of 
molten  KHSO4  (See  reference  [49]  and  references  therdn) 

accotmted  for  in  the  calculations.  However,  the  qualitative  agreement  is  good,  the 
pattern  of  vibrational  levels  bdng  well  reproduced.  For  the  intensities  there  are  no 
absolute  measurements.  However,  it  is  known  that  1^  fr  the  most  intense  IR  band  for 
SO^**,  and  that  ui  is  the  most  intense  Raman  band,  which  is  dearly  demonstrated  in 
the  theoretical  resiilts. 

The  computed  harmonic  frequendes  and  intensities  of  HSOJ  are  p ven  in  Thble  2. 
We  denote  them  according  to  thdr  symmetry  label  in  the  C,  point  group.  We  also 
present  experimental  frequendes  for  HSO4  in  the  NH4HSO4  crystal,  taken  frx>m  the 
work  of  Dawson  et  al.  [49].  We  see  overall  quatitative  agreement  between  the  theoret¬ 
ical  and  experimental  frequendes,  with  the  S-0  stretdies  bdng  somewhat  too  high 
as  found  for  SOU',  except  for  the  S-OH  strdch  wlddi  is  too  low,  possibly  correlating 
with  the  error  found  for  the  S-OH  bond  distance,  which  was  too  long.  Quantitative 
agreement  should  not  be  expected  due  to  both  limitations  in  the  theory  and  other 
effects  such  as  additional  hydrogen  bonding  in  the  solid,  which  is  expecAed  to  lower 
the  experimental  0-H  stretching  frequency.  As  far  as  the  intensities  are  conconed, 
the  computed  results  seem  to  be  in  fur  agreement  vnth  the  experim€aitally  observed 
values  in  that  the  most  intense  Raman  band  experimentally  is  the  1022-1044  cm'* 
band  and  this  corresponds  to  the  1157  cm'*  band  in  the  theoretical  calculations. 
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Figure  3:  structure  in  the  presence  of  an  electric  field  along  a  C3  axis.  Bond 

lengths  are  in  a.u.  and  angles  are  in  degrees  (for  the  electric  field,  1  a.u.  = 
5.14221x10®  V/cm). 

Fhothermore  the  order  and  relative  intensity  of  the  four  bands  at  642,  862,  1157, 
and  1258  cm*'^  axe  not  expected  to  change  because  aU  these  modes  have  the  same 
symmetry.  Overall,  the  theoretical  methods  used  here  give  useful  qualiUdiive,  though 
not  quantitative,  results. 

4.2  S04~  in  the  presence  of  an  electric  field 

In  this  section  we  consider  the  effect  of  an  implied  electric  field  on  the  structure, 
frequencies  and  IR  intensities  of  SO®'.  This  analyris  is  useful  as  a  simple  model  of 
electric  field  effects  at  an  electrode  surface,  and  also  in  the  analysis  of  binding  in  the 
copper  plus  sulfate  clusters.  The  magnitude  of  the  field  was  fixed  at  0.01  a.u.  (5x10®^ 
V/cm),  corresponding  to  the  strong  electric  field  expected  in  the  electrochemical 
double  layer  (13, 50]. 

Ihe  effect  of  the  iq>plied  fidd  on  the  structure  of  SO®'  is  illustrated  in  Figure  3. 
The  field  is  applied  along  a  Cz  asds,  with  a  poritive  and  n^ative  rign,  and  the  SO®' 
unit  has  Cs«  symmetry  in  b<^  cases,  with  different  bond  distances  and  axigles.  Con¬ 
sidering  first  Figure  3(a),  we  see  that  the  S-*01  bond  has  lengthened  significantly  frmn 
the  isolated  ion  value  of  2.80  a.u.,  whereas  the  S-02  bond  has  shortened.  However, 
the  shortening  effect  is  smaller  rince  the  S-02  bond  lies  at  a  greater  angle  to  the  iq>- 
plied  field.  The  bond  an^es  show  effects  consistent  with  the  distances.  The  01-S-02 
angle  has  decreased  from  the  iadated  ion  value  of  109.5**,  whereas  the  02-S-02  angle 
has  increased.  Thinking  of  the  S  position  as  fixed,  we  see  that  all  the  O  atoms  have 
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Tabic  3:  Summary  of  the  SO^  frequencies  (cm"')  and  IR  intensities  (km/mol)  with 
an  electric  field  along  S-Ol  (Cs^  symmetry) 


-0.01*  a.u. 

40.01*  a.u.  1 

V 

1  IRint. 

u 

IRint. 

l/2(c) 

471 

0.1 

471 

0.1 

666 

54.2 

665 

44.6 

664 

37.8 

667 

58.8 

*'i(cii) 

1015 

138.2 

1028 

20.3 

l/3(c) 

1210 

583.0 

1136 

626.3 

*^(ai) 

1113 

517.3 

1248 

547.5 

*  The  orientation  of  the  field  corresponds  to  Figure  3(a) 
^  The  orientation  of  the  field  coiresponds  to  Figure  3(b) 


moved  in  the  direction  opposite  to  the  applied  field,  as  may  be  expected  due  to  thdr 
negative  charge.  Reversing  the  direction  of  the  applied  field  (Figure  3(b))  has  the 
expected  effect  —  the  O  atoms  all  move  in  the  opposite  direction,  so  that  the  bond 
distances  and  angles  change  in  the  opporite  sense  but  by  about  the  same  magnitude 
as  for  the  first  field  direction. 

The  computed  harmonic  frequencies  and  IR  intensities  for  in  «he  presence 
of  an  applied  electric  field  are  given  in  Thble  3,  which  are  to  be  compared  with  the 
isolated  ion  results  of  Table  1.  The  frequencies  in  Table  3  are  identified  by  the  original 
labels  (t/j  —  1/4)  used  in  Table  1,  as  the  assignment  is  quite  clear  even  though  the 
symmetry  is  lower.  We  see  that  the  1%  modes  of  the  Ti  point  group  have  split  into  ai 
and  e  modes  for  the  Czv  point  group  in  the  presence  of  the  applied  field.  Also,  uz  and 
V4  are  barely  affected  by  the  applied  field,  presumably  because  these  modes  involve 
motions  which  are  predominantly  perpendicular  to  the  applied  fidd.  Because  of  their 
low  IR  and  Raman  intensity,  these  modes  are  difficult  to  detect  experimentally. 

Of  more  interest  are  vi  and  as  these  are  intense  in  the  Raman  and  IR  r^ons, 
respectively  (see  Table  1).  Although  we  have  not  computed  Raman  intensities  in 
the  presence  of  an  applied  field,  it  is  likely  that  ui  will  still  be  the  most  int^ise 
Raman  band  in  this  case  also.  From  Table  3,  for  the  -0.01  a.u.  fidd,  we  see  that  the 
lengthening  of  the  S-Ol  bond  (Figure  3(a))  results  in  a  reduction  of  the  i^(ai)  mode 
by  about  20  cm**'  and  the  f^(oi)  mode  by  about  60  cm"',  a  large  effect.  In  addition, 
Vi  has  gained  some  IR  intensity,  possibly  due  to  interaction  with  the  intense  i^(ai) 
mode.  The  t^(e)  mode  has  increased  in  frequency  by  about  35  cm"',  reflecting  the 
shorter  S-02  bond  distance.  We  note  that  pi  involves  both  S-Ol  and  S-02  stretches, 
so  that  in  the  presence  of  the  applied  field  there  is  a  balance  between  effects  which 
increase  u\  (shortening  the  three  S-02  distances)  and  decrease  i/j  (lengthening  the 
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S-Ol  distance). 

From  the  +0.01  a.u.  results  of  Table  3,  we  sec  that  reversing  the  direction  of 
the  applied  field  has  a  similar  effect  on  the  frequencies  as  for  the  bond  distances. 
The  Ui(ai)  mode  is  again  reduced  in  frequency,  although  the  effect  is  smaller  than 
for  the  -0.01  a.n.  field.  This  shows,  however,  that  an  intuitive  guess  based  on  the 
bond  distances  is  not  entirely  appropriate  for  predicting  the  frequency  shifts  —  in 
this  case,  the  additional  us{ai)  mode  has  an  effect  on  the  i^(ai)  mode,  in  one  case 
moving  down  in  frequency  and  the  other  case  moving  up.  It  is  probable  that  this 
shift  in  i^(aj)  has  a  subtle  but  important  effect  on  for  the  -0.01  a.u.  field  case, 
increasing  the  IR  intensity  and  reducing  the  frequency.  FWther  evidence  for  this  may 
be  seen  from  the  IR  intensity  of  in  the  +0.01  a.u.  field  case  —  the  increase  is 

much  smaller  than  for  the  -0.01  a.u.  case,  because  i<s(ai)  has  moved  up  in  frequency 
rather  than  down.  The  *^(e)  modes  are  shifted  down  in  frequency  for  the  +0.01  a.u. 
field  case,  as  expected. 

If  we  were  to  consider  the  metal  surface  to  be  perpendicular  to  the  applied  field 
direction,  then  we  must  keep  in  mind  the  fact  that  the  c  modes  will  not  be  observable 
on  a  surface.  This  is  because  the  IR  intensity  arises  from  a  dipole  moment  derivative 
which  is  parallel  to  the  “surface”,  which  will  be  cancelled  by  image  charge  effects  (ie. 
the  surface  selection  rule).  Therefore  it  is  the  ai  modes  which  are  the  most  important, 
as  these  involve  dipole  moment  derivatives  which  ate  perpendicular  to  the  “surface” . 
These  simple  model  calculations  indicate  that  the  orientation  of  the  S04~  molecule 
on  the  surface  may  be  determined  by  the  fr^sqnen^  shift  of  the  intense  i^(ai)  mode. 

4.3  Cu^SOj"  results 

In  this  section  we  consider  two  Cu^SO’”  structures,  illustrated  in  Figure  4,  denoting 
them  as  *l-down”  and  “3-down”,  respectively.  There  are  other  structures  which  may 
be  important  —  for  example,  there  is  a  *1.5  down"  structure  with  a  Cu-Ol-S  angle 
of  122*.  This  lies  around  2  kcal/mol  bdow  the  1-down  structure.  In  the  current  work 
we  consider  only  the  Cjv  structures,  postponing  discussion  of  alternative  structures 
to  future  work  [19] 

The  interaction  of  SO*"  with  the  Cu  atom  is  mainly  electrostatic,  dominated 
by  charge-induced  dipole  interactions.  Remembering  that  SO*  has  no  permanent 
dipole  moment,  this  means  that  the  binding  energy  of  SO*"  to  a  single  Cu  atom  is 
not  especially  sensitive  to  the  orientation  of  the  SO*"  unit.  The  4s  electron  on  Cu 
polarizes  away  from  the  SO*"  unit,  forming  an  sp  hybrid  orbital  in  order  to  both 
reduce  the  repulsion  and  create  an  attractive  riectrostatic  interaction  by  forming  a 
positive  re^on  of  charge  on  the  SO*"  side  of  the  atom  —  an  induced  dipole.  The 
interaction  is  strong  —  the  binding  energy  of  the  structure  illustrated  in  Figure  4(a) 
is  around  37  kcal/mol.  For  comparison,  at  an  equivalent  level  of  theory  the  binding 
energy  of  Cu°H20  is  less  than  1  kcal/mol  (a  Van  dcr  Waals  complex),  whereas  the 
binding  energy  of  Cu'^HjO  is  strong  at  around  32  kcal/mol. 
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S-02  2.77  Z.Cu°-02-S  79* 


(a)  (b) 


Cs*  (*1  down*);  (’3  down*); 

2  kcal/mol  8  kcol/mol 

Figure  4:  Summary  of  Cu^’SO^'  structures  in  Ca«  symmetry.  Bond  lengths  axe  in 
a.u.  and  angles  axe  in  degrees. 

The  two  structures  are  separated  energetically  by  6  kcal/mol.  The  Cu>0  bond 
distances  in  each  case  are  very  different.  The  Cu-Ol  distance  the  1-down  structure 
is  reasonable  —  for  example,  the  eariy  surface  exa£i  experiments  (5]  found  a  value  of 
3.9  a.u.  The  0u-O2  distance  of  4.83  a.u.  for  the  3-down  structure  is  much  greater  — 
in  this  case  the  limiting  factor  is  the  inal^ty  of  the  S04~  unit  to  distort  sufficiently 
to  allow  a  stronger  interaction  between  the  02  atoms  and  the  Cu,  combined  with  the 
repulrion  between  the  formally  poritivriy  charged  S  and  the  poritive  charge  due  to 
the  induced  dipole  on  the  Cu  atom.  Hie  fact  that  there  axe  three  02  atoms  involved 
in  the  interaction  offsets  these  factors  to  some  extent. 

It  is  of  interest  to  compare  the  bond  lengths  and  angles  of  the  SO^"  unit  in 
Cu^O^**  (Figure  4)  with  those  of  S04~  in  the  presence  of  an  ^plied  firid  (Figure  3). 
Fbr  the  3-dowtt  structure,  we  note  that  the  S-02  distance  is  the  same  as  in  free  SO^", 
in  contrast  to  Figure  3(b),  whoe  it  is  slightly  longer.  However,  the  02-S-O2  angle 
is  smaller  for  the  3-down  structure  than  in  Figure  3(b).  These  differences  are  due 
to  the  more  localized  nature  of  the  electrostatic  interactkm  betweoi  Cu  and  SOU* 
ctnnpared  vdth  the  uniform  electric  field  of  Figure  3.  Overall,  however,  tiMie  is  a 
remarkable  similarity  between  the  SOl"  unit  in  the  1-down  structure  and  the  applied 
field  structure  of  Figure  3(a),  and  the  SO4'  unit  in  the  3-down  structure  and  the 
applied  field  structure  of  Figure  3(b).  Tins  is  due  to  the  dominantly  electrostatic 
nature  of  the  interaction  of  S04~  with  Cu  —  the  Cu  atom  with  an  induced  dipole 
is  quite  similar  to  an  applied  electric  firid,  as  far  as  the  SOl"  unit  is  amconed.  In 
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S-01  2.73 

S-02  3.04  S-02  3.03 

S-03  2.7t  S-03  2.70 

Z.03-S-03  114.6*  /.03-S-03  114.6* 


(a)  (b) 

Qi  ('1  down*);  C»  (*1  down*); 

Okccl/mot  <0.5  kcal/mol 

Figure  5:  Summaiy  of  Cu^HSOj  structures  in  C,  symmetry.  Bond  lengths  are  in 
a.u.  and  angles  are  in  degrees. 

future  work  [19]  we  will  consider  whether  this  also  implies  to  the  frequenqr  shifts 
for  the  SO^”  unit  in  Cu’^SO^**.  At  the  present  time  we  amply  state  that  they  are 
expected  to  be  qualitatively  similar  to  those  found  for  S04~  in  the  presence  of  an 
applied  field. 

4.4  Cu®HS04  structure 

We  now  consider  two  Cu^HSOj  structures,  illustrated  in  Figure  5.  As  for  Cu^SO^'", 
there  are  other  possible  structures.  However,  these  are  the  lowest  energy  structures 
which  we  have  found. 

Hiere  are  two  important  differences  between  and  HSO4 ,  both  due  to  the 
presence  of  the  H  in  HSOJ.  The  first  is  the  reduced  diarge  of  HSOJ,  whidbi  lesuits  in 
a  smaller  diarge-induoed  dipole  interaction  in  Cu^HSOJ.  This  is  illustrated  by  the 
reduced  binding  enogy  of  Cu^HSO^v  only  around  14  kcal/mol,  less  than  half  that 
of  Cu**S04~.  The  second  difference  is  the  presence  of  a  permanent  dipole  moment 
in  HSO7.  Considoing  Figure  2(b),  the  permanent  dipole  is  oriented  in  the  S>02-H 
plane,  approximately  perpendicular  to  the  03-S-03  plane,  pasang  through  the  S 
with  the  positive  end  at  the  OH  end  of  the  molecule.  There  will  be  an  effect  which 
tends  to  orient  the  dipole  to  maxinuxe  the  dipole-induced  dipole  interaction.  However, 
this  is  a  secondary  dfect  compared  with  the  charge-induced  dipole  interaction.  This 
is  illustrated  by  the  two  structures  shown  here  —  the  orientation  of  the  dipole  is 
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different  in  each  case,  yet  the  structures  have  very  similar  energies. 

The  Cu-Ol  bond  distance  is  longer  in  Cu^HSO^  than  Cu°S04~,  reflecting  the 
lower  binding  energy  of  Cu^HSOi" .  Comparing  the  HSO4  structure  in  Cu®HSO^ 
with  the  isolated  HSO7  structure  of  Figure  2(b),  we  see  similar  effects  to  the  SO^” 
structural  changes  in  Cu^SOl",  but  again  these  are  not  as  large  since  the  interaction 
is  weaker.  In  addition,  the  symmetry  is  lower.  The  S-Ol  bond  distance  is  greater 
in  Cu^HSO^  compared  with  HSO^*  and  the  S-02  and  S-03  distances  are  shorter. 
The  angles  also  change  in  a  way  consistent  with  that  found  for  SO^'  in  Cu®S04~  — 
the  negatively  charged  oxygens  move  toward  the  positive  dde  of  the  Cu,  whereas  the 
positively  charged  H  moves  away.  This  decreases  the  Ol-S-02  angle,  for  example, 
and  increases  the  S-02-H  angle. 

Th^  are  several  other  possible  structures  for  Cu®HS04  when  compared  to 
Cu®SOj“,  mainly  because  of  the  lower  symmetry  of  HSO4  compared  to  SOj".  How* 
ever,  the  most  important  difference  between  Cu®S04~  and  Cu®HS07t  that  is  the 
weaker  binding  leading  to  smaller  geometric  distortions  of  HSO7  in  Cu^HSO^ ,  will 
also  apply  in  these  other  cases.  This  also  means  that,  if  HSO7  and  SOl"  adsorb  on  a 
Cu  surface  with  a  similar  orientation,  the  frequency  shifts  for  HSO7  may  be  expected 
to  be  smaller  than  for  SOl".  However,  since  the  frequency  shifts  are  very  dependent 
on  the  orientation,  this  cannot  be  taken  as  a  general  rule. 

4.5  CuJSOj”  structures 

In  this  section  we  present  results  for  two  structures  of  Cu^SO^',  illustrated  in  Figure  6. 
These  are  chosen  to  illustrate  important  liimtations  in  the  single  copper  atom  model, 
but  also  to  illustrate  the  utility  of  that  model  in  the  qualitative  understanding  of  the 
interaction  of  S04~  and  a  copper  surface. 

One  valid  question  is  why  consider  the  single  copper  atom  calculations  at  all, 
given  that  we  are  also  presenting  the  Cuf  results.  There  are  several  reasons  —  firstly, 
it  is  computationally  fearible  to  compute  harmonic  frequencies  and  IR  intensities  for 
Cu^SO’"  and  related  systems.  Although  we  have  not  presented  these  results  here, 
this  will  appear  in  future  work  [19].  Secondly,  the  Cu®S04~  results  give  good  insight 
into  the  nature  of  the  CU-SO4''  interaction,  without  the  additional  complications  of 
duster  artifacts  which  may  ^pear  (see  below).  Finally,  the  single  copper  calculations 
serve  as  a  starting  point  from  which  to  plan  the  more  computationally  intensive  CU4 
calculations. 

The  four  atom  copper  duster  used  here  illustrates  the  additional  binding  possi¬ 
bilities  of  S04~  on  the  Cu(lll)  surface  —  the  3^own  and  1-down  models  now  contain 
interactions  with  ihrtie  copper  atoms.  Before  discussing  Cu^SO^",  however,  we  must 
note  some  ad<iition.al  limitations  of  the  duster  modd. 

For  the  CU4  and  the  Cu^SO^**  dusters,  we  have  currently  used  a  dosed-shell 
singlet  coupling  of  the  Gu  4s  electrons.  There  are  other  possible  couplings  —  in 
particular,  several  triplet  states,  which  we  will  be  studying  in  future  work.  Whether 
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Figure  6:  Summary  of  Cu^SOl*  structures  showing  ontop  and  three  fold  hollow 
sites.  Bond  lengths  are  in  a.u.  and  angles  are  in  d^ees. 


a  nnglet  or  triplet  coupling  the  electrons  is  duMen,  there  are  problems  which  arise 
from  unng  the  cluster  model.  Cu4  is  in  reality  a  small  molecule  with  a  valenoe  shell  of 
4s>like  electrons,  rather  than  an  accurate  model  of  a  metal  surface.  The  occupation  of 
the  valence  shell  leads  to  an  asymmetric  charge  distribution  in  symmetry,  so  that 

Cul  in  Figure  6  carries  a  net  negative  diarge  and  the  three  Cu2  atoms  carry  an  equal 
net  positive  charge,  for  the  closed-shell  singlet  coupled  case.  This  gives  the  cluster 
an  overall  dipole  moment  which  in  this  orientation  favours  the  binding  of  a  n^ative 
ion  such  as  SOI',  and  so  absolute  binding  ener^es  are  not  well  defined.  However, 
relative  binding  ener^es  should  still  be  qualitatively  useful. 

The  two  structures  of  Figure  6  illustrate  the  importance  of  multidentate  bond¬ 
ing  for  SO|'  on  copper.  The  S-down  structure  is  mu^  more  stable  than  the  l-down 
structure  over  the  3-fold  hollow  ate.  The  binding  energur  of  S0|'  to  Cu*  (dosed  shell 
singlet  coupled)  for  the  3-down  structure  is  54  kcal/mol,  compared  with  31  kcal/mol 
for  the  1-down  stru<^ure.  The  binding  energy  for  the  3-down  structure  is  also  much 
hi^er  than  that  of  S0|'  bound  to  Cu  (37  kcal/mol,  §  4.3).  Some  part  of  the  increased 
binding  energy  for  the  Cu4  duster  in  the  3-down  orientation  is  due  to  the  permanent 
dipde  moment  of  the  Cu<  duster.  However,  the  large  difference  of  23  kcal/mol  be¬ 
tween  the  3-down  and  1-down  cases  for  Cu|SO|'  indicates  that  the  3-down  structure 
may  be  favoured  on  Cu(lll).  Interestingly,  the  Cu-0  distances  are  very  similar  for 
both  cases,  just  over  4.1  a.u.  However,  the  interaction  of  three  n^ativdy  charged 
O  atoms  in  the  3-down  case  obviously  leads  to  a  much  stronger  interaction  than  the 
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l-<iown  case.  The  Cu-0  bond  ien^h  of  4.1  a.u.  is  reasonable  when  compared  with 
the  surface  exafs  value  of  3.9  a.u.  discussed  previously.  However,  we  do  not  expect 
quantitative  agreement. 

It  is  useful  to  compare  the  structures  of  the  units  in  Figure  6  with  the  ear¬ 
lier  Figures  3  and  4.  The  SO^”  unit  for  the  3-down  structure  of  Figure  6(a)  is  quite 
similar  to  that  of  Figure  3(b),  indicating  again  a  dominantly  electrostatic  interaction. 
Compared  to  the  Cu^SOj"  result  (Figure  4(b)),  the  angles  of  the  SOj"  unit  in  the 
Cu4S04~  3-<lown  structure  indicate  less  distortion,  due  to  the  oxygens  being  in  a  more 
favourable  position  in  the  CujSO*"  structure.  The  l-down  CujSOj*"  structure,  al¬ 
though  mu^  higher  in  energy,  is  still  interesting.  We  see  again  qualitative  similarities 
to  the  earlier  results  of  Figures  3(b)  and  4(b).  However,  the  distortions  of  the  ari¬ 
ses  and  the  S-Ol  bond  distance  are  mudi  greater  in  the  l-down  Cu^SOl'structure. 
Although  the  interaction  is  probably  still  dominantly  electrostatic  for  the  1-down 
structure,  the  very  long  S-Ol  bond  distance  may  be  indicative  of  the  formation  of 
a  weak  chemical  ^nd  between  the  Gui  duster  and  the  unit.  Comparing  the 
bond  lengths  and  angles  of  the  CuJSO*”  1-down  structure  with  those  of  EIS07  (Fig¬ 
ure  2(b)),  we  see  qualitative  similarities.  Althought  the  distortions  are  not  as  large 
in  Cu^SO^^  as  in  HSOJ  it  is  dear  that  the  1-down  structure  is  some  way  along  to  a 
chemically  bonded  spedes,  in  contrast  to  the  3-down  structure,  which  is  dectrostati- 
cally  bound. 

Regarding  fr«qucn<y  shifts  for  SO’"  in  CuJSO’",  we  simply  note  from  the  geo¬ 
metric  distortions  of  the  SO’"  units  that  they  axe  expected  to  be  qualitativdy  similar 
to  those  found  for  S0’~  in  the  presence  of  an  ^plied  fidd,  and  the  1-down  shifts  may 
be  larger  than  the  3-down  shifts. 


5  Conclusions 

A  series  of  modd  calculations  have  been  presented  which  aim  at  hdping  to  interpret 
the  results  of  in  situ  experimental  data  on  the  adsorption  of  SO’"  and  HSO4  on 
metal  surfaces.  The  current  calculations  are  aimed  at  adsorption  on  copper,  but  the 
qualitative  nature  of  the  results  means  that  the  condusions  may  generaUze  to  otho: 
metals.  The  geometry  optimization,  harmonic  frequency  and  intensity  calculations 
for  isolated  SO’"  and  HSO4  indicate  useful  qualitative  accuraQr.  The  application  of 
a  strong  uniform  dectric  field  to  SO’"  gave  insight  into  the  response  of  SO’"  to  the 
local  field  in  the  dectrodiemical  double  layer,  and  also  gave  insight  into  the  nature 
of  the  bonding  of  SO’"  to  Cu  and  OU4.  The  Crequenqr  shifts  computed  for  SO’"  in 
the  presence  of  the  applied  field  could  be  straightforwardly  interpreted  in  terms  of 
the  geometrical  distortions  caused  by  the  fidd.  However,  the  calculations  show  that 
there  are  some  subtle  effects  for  the  frequency  shifts  due  to  the  interaction  of  modes 
of  the  same  symmetry  which  are  not  obvious  from  intuitive  arguments  based  on  the 
geometrical  distortions.  In  addition,  the  calculations  show  that  the  frequency  shift  of 
the  most  intense  «!  mode  may  serve  to  indicate  the  adsorption  geometry  of  SO’"  on 
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Cu(lll),  assiuning  symmetry. 

The  interaction  between  Cu®  and  SOj'  or  HSO4  is  dominantly  ionic,  a  charge- 
induced  dipole  interaction,  so  there  is  not  a  large  energetic  difference  between  the 
different  structures.  The  binding  energy  of  SO*”  to  Cu  or  Cu*  is  quite  high,  whereas 
the  binding  energy  HSO7  to  Cu  is  lower  due  to  the  reduced  charge  of  HSO^i  giving  a 
reduced  charge-induced  dipole  interaction.  The  similarity  of  the  structure  of  the  SO^” 
unit  in  the  presence  of  an  iqiplied  field  with  that  in  Cu®S04”  serves  to  illustrate  the 
ionic  nature  of  the  bonding.  The  lower  binding  energy  of  HSO7  compared  with  SO*” 
means,  that  for  an  equivalent  adsorption  orientation,  the  frequency  shifts  observed 
for  HSO4  may  be  less  than  for  SO*”.  However,  since  the  frequency  shifts  are  very 
dependent  on  adsorption  geometry,  this  statement  cannot  be  taken  as  general. 

The  Cu®SO*~  calculations  illustrate  that,  ev^  accounting  for  artifacts  due  to 
the  cluster  model,  the  3-down  multidentate  structure  of  Cu^SO*”  is  energetically 
much  more  favourable  than  the  1-down  hollow  rite  adsorption  for  Cu(lll).  This  is 
mainly  due  to  the  very  favourable  position  of  the  oxygen  atoms  over  the  Cu  atoms 
on  Cu(lll).  The  comparison  of  the  SO*”  structure  in  CujSO*”  with  the  earlier 
structures  indicates  that  the  qualitative  nature  of  the  frequency  shifts  for  SO^”  may 
be  deduced  from  the  geometric  distortions  of  the  SO*”  units.  Finally,  the  structure 
of  the  SO*”  unit  in  the  1-down  CujSO*”  cluster  indicates  the  possible  formation  of 
a  weak  chemical  bond  when  compared  to  the  3-down  Cti^SO*”  structure,  whidi  is 
electrostatically  bound. 

Continuing  and  future  work  includes  the  consideration  of  more  structures  for 
Cu®SO*”  and  Cu^HSO^ ,  the  computation  of  harmonic  frequendes  and  IR  intensities 
for  these  spedes,  the  optimization  of  additional  structures  for  the  larger  Cu^SO*" 
structures  and  the  consideration  of  even  larger  dusters  such  as  Cu®oSO*~.  When 
combined  with  the  current  results,  these  should  give  more  insight  into  the  interaction 
of  SO*”  and  HSO7  with  metals  and  metal  surfaces. 
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